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Abstract 

Quantum ergodicity, which expresses the semiclassical convergence of almost all 
expectation values of observables in eigenstates of the quantum Hamiltonian to the 
corresponding classical microcanonical average, is proven for non-relativistic quantum 
particles with spin 1/2. It is shown that quantum ergodicity holds, if a suitable 
combination of the classical translational dynamics and the spin dynamics along the 
trajectories of the translational motion is ergodic. 
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1 Introduction 



In quantum chaos one is primarily interested in statistical properties of eigenvalues and 
eigenvectors of quantum Hamiltonians whose classical limits generate chaotic dynamics. 
In this context the eigenvalue statistics on the scale of the mean level spacing are expected 
to be described by random matrix theory (see, e.g., [|BGS84| , |Meh91|] ). This conjecture 
has since found overwhelming confirmation, mostly based on numerical calculations of 
eigenvalues in many systems. For the corresponding eigenvectors one also expects a random 
behaviour, and various tools to measure this have been invented. Among them one finds one 
of the few results in quantum chaos for which a mathematical proof is available: if a classical 
system is ergodic, its quantum mechanical counterpart is quantum ergodic. By this one 
understands a semiclassical convergence of almost all phase-space lifts of eigenfunctions of 
the Hamiltonian (e.g., their Wigner- and Husimitransforms) towards equidistribution with 
respect to microcanonical (i.e. Liouville) measure. One thus has obtained a realisation of 
the semiclassical eigenfunction hypothesis (see, e.g., ||Ber83|| ) for classically ergodic systems. 

Quantum ergodicity was first established for the free motion of a particle on a compact 
Riemannian manifold, where the quantum dynamics is generated by minus the Laplace- 
Beltrami operator on that manifold. This goes back to Shnirelman ||Shn74|] , and the first 
complete proofs are due to Zelditch [|Zel87|| and Colin de Verdiere ||CdV85|| . In the systems 
considered by these authors the semiclassical limit is actually realised as a high-energy limit. 
However, for Schrodinger operators involving a potential and possibly also a magnetic field 
the semiclassical limit can in general only be performed in terms of h 0. In this setting 
quantum ergodicity was first proven by Helffer, Martinez, and Robert |[IIMR87| . All of 
the work mentioned so far exclusively dealt with the dynamics of point particles without 
internal degrees of freedom; in particular no spin was involved. On the level of the quantum 
Hamiltonians H that were taken into account, this is reflected in the fact that H appears 
as a (Weyl-) quantisation of a scalar symbol. The only degrees of freedom appearing are 
therefore those that possess a direct classical analogue. These are in fact the translational 
degrees of freedom, which on the classical side yield as a phase space a smooth symplectic 
manifold that in typical cases is the cotangent bundle over the configuration manifold. In 
this setting internal degrees of freedom would appear through the fact that the quantum 
mechanical observables are quantisations of matrix valued symbols, such that the vector 
spaces these matrices operate on represent the internal degrees of freedom. 

In quantum mechanics internal degrees of freedom often arise due to the presence of 
symmetries. Since in general symmetries have to be implemented through (anti-) unitary 
representations of the respective symmetry groups on the Hilbert space of state vectors, 
representations of dimensions exceeding one introduce discrete degrees of freedom. It then 
can happen that these additional degrees of freedom do not possess classical analogues. A 
prominent example of this phenomenon is provided by the spin of a particle, which arises 
through the space-time symmetries that either form the Lorentz group (in a relativistic 
theory), or the Galilei group (in a non-relativistic theory). Both these symmetry groups 
contain the proper rotations, i.e. SO (3), as a subgroup. Through the unitary projective 
representations of the Lorentz or Galilei group, respectively, which implement the space- 
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time symmetries in the quantum theory, one therefore introduces the unitary irreducible 
representations of Spin (3) = SU(2), i.e. spin. On the classical side the rotation group, 
however, remains to be SO (3), which then enters in terms of spatial angular momentum. 
There exists, however, no direct classical analogue of spin. Coming back now to the problem 
of quantum ergodicity for particles with spin, there immediately arises the question for a 
criterion to be imposed on the classical system in order to ensure that the eigenfunctions 
of the quantum Hamiltonian behave quantum ergodically in the semiclassical limit. We 
recall that in the case without spin such a criterion was given by the ergodicity of the 
dynamics on the classical phase space with respect to Liouville measure. One would now 
expect that the presence of additional, internal degrees of freedom requires an extended 
criterion in order that quantum ergodicity holds. 

It is the primary goal of the present work to elaborate on this question in some detail 
for the case of a non-relativistic quantum particle with spin 1/2. Our proof of quantum 
ergodicity in this context generalises the methods of ||Shn74| , |Zel87| , |(JdV85| , [H1V1K87[ to the 
situation of Weyl operators with 2 x 2-matrix valued symbols. This requires two essential 
ingredients. The first one is an Egorov Theorem, which relates the semiclassical limit of 
the quantum mechanical time evolution of an observable to the classical time evolution of 
the corresponding classical observable. The second input required is a Szego limit formula, 
which expresses the semiclassical limit of averaged expectation values of observables in 
eigenstates of the Hamiltonian in terms of a classical microcanonical average. In the course 
of the subsequent proof of quantum ergodicity we primarily rely on the method developed in 
Zcl96| , |ZZ96|| , which provides a considerable simplification over the original proofs given in 



Zel87| , |CdV85| , [HMK87|| . Already the Egorov property that one encounters for systems with 



spin hints at the construction that yields the 'classical' criterion for quantum ergodicity 
that we seek for. It leads us to consider an SU(2)-extension of the Hamiltonian flow that 
arises from the classical limit of the translational degrees of freedom. The thus extended 
flow is defined on a product phase space that consists of two parts: the translational part 
is the hypersurface of fixed energy in the classical phase space, whereas the spin part is 
given by the group manifold of SU(2). The translational part of the combined dynamics 
is then provided by the Hamiltonian flow. The latter also drives the spin dynamics that 
takes place on SU(2) and consists of a left multiplication by a spin transport matrix, which 
propagates the spin degrees of freedom along the trajectories of the Hamiltonian flow. 
Our main result, stated in Theorem [4.1| , then is that ergodicity of the combined flow is 
a sufficient criterion for quantum ergodicity. At this point we remark that the problem 
of quantum ergodicity for quantisations of matrix valued symbols was already considered 
Zel96|| as an example for the general theory of quantum ergodicity of C* dynamical 
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systems. There it was, however, overlooked that in the relevant Egorov Theorem, apart 
from a transport by the Hamiltonian flow, the principal symbol of an observable is also 
conjugated with the spin transport matrices. As a consequence, the result stated in ||Zel96 | 
is hence incorrect, to the extent that ergodicity of the Hamiltonian flow alone is insufficient 
to guarantee quantum ergodicity. In section ^we give an example that illustrates this fact. 

This paper is organised as follows. In section ^ we review some background on h- 
pseudo differential calculus and define the type of Hamiltonians and observables that will 
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be considered in the sequel. The two major ingredients required for quantum ergodicity, 
namely the Egorov property and the Szego limit formula, are developed in section ^. Our 
main result, quantum ergodicity for Pauli Hamiltonians with spin 1/2, is then proven in 
section |. Finally, in section ^ we discuss the consequences of quantum ergodicity for 
Wigner- and Husimitransforms of eigenfunctions and, furthermore, give an example that 
illustrates why ergodicity of the Hamiltonian flow alone is not a sufficient criterion for 
quantum ergodicity. 



2 Semiclassical background 

In non-relativistic quantum mechanics the dynamics of a particle with spin s G |N is 
governed by the Pauli equation 

ih^{x,t) = Hp^{x,t) (2.1) 

with the quantum Hamiltonian 

Hp = HtransX2s+l + ■ C (2.2) 

acting as a self-adjoint operator on a suitable domain in the Hilbert space L^(]R'^) ® C^*+^. 
Here 

I fh e V 
Htrans. = — "V^ - -A{x) + e^{x) (2.3) 
2m \i c J 

describes the dynamics of the translational degrees of freedom of a spinning particle with 
mass m and charge e which is subject to external electromagnetic forces generated by the 
(static) potentials A and (p. Furthermore, the components S^, k = 1,2, 3, of S denote 
the hermitian generators of the Lie algebra su(2) in the (irreducible) spin-s representation, 
which is of dimension 2s+l. The coupling of the spin degrees of freedom to the translational 
ones is provided by the operators Ck, which are suitable quantisations of functions Ck{p, x) 
on phase space. The latter can, e.g., describe a coupling to an external magnetic fleld, 
Cb{p,x) = -^^B{x), or a spin-orbit coupling Cso{p,x) = ^^■^l^^^i'^'fff {x x p). 

Here, and in the following, we choose all quantum mechanical observables to be Weyl 
quantisations of matrix valued symbols. In general this is deflned for B G iS'(M'^ x M'^) 
C"^" and G 5(M'^) ® C" as 

(5^)(x) := -i— / / e^^<'^-y'>Bfp,]-{x + y))^{y)dydp, (2.4) 
[Ann) j^d j^d \ A / 



where h G {0, Hq] serves as a parameter. The quantisation yields a continuous map 
from 5(]R^) Cg) C" to i5'(]R^) C". However, in order to obtain a semiclassical calculus one 
has to restrict attention to smaller classes of symbols and operators. In particular, one 
wishes to consider operators that can be composed with one another, e.g., operators which 
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map 5(M'') ® C" into itself. In defining suitable symbol classes we use the following notion, 
which is in accordance with [PS99|| : 



A function m : M*^ x M'^ ^ (0, oo) is called an order function, if there are constants 
Co > 0, No>0 such that 

m{p,x) < Cqa/I + {p-qY + {x- yf^" rn{q,y) . (2.5) 
An example for such an order function is 

m{p, x) = 1 + p"^ + x'^ . (2.6) 

Definition 2.1. Let m be an order function on x W^. We define the symbol class S(m) 
to be the set of S G C°°(R"' x M'^) ^C"""" such that for every a, /? G N[J there exists Ca,f3 > 
with 

\\d^d^,B{p,x)\\<C^,pm{p,x) , (2.7) 
where || ■ || is some (arbitrary) matrix norm on C"^". 

If i? = B{p, x; h) depends on h E (0, 1], we say that B G S(m), if i?(-, ■; h) is uniformly 
bounded in S(m) when h varies in (0, 1]. For /c G M we let S'^(m) be the set of functions 
B{jp, x; h) on M'^ x M'^ x (0, 1] that belong to h~^S{m) and satisfy 

\\d;dlB{p,x-h)\\ < C^,^m{p,x)h-' . (2.8) 

A sequence of symbols Bj G S^j (m), with kj — oo monotonically, defines an asymptotic 
expansion of i? G S'^"(m), denoted by 

oo 

, (2.9) 

j=0 

if 

N 

B -J^Bj eS''''+'{m) (2.10) 

j=0 

for every G Nq. We will often use the smaller class of classical symbols S^i(m), whose 
elements B G S^j(m) possess asymptotic expansions in integer powers of h, i.e. 

oo 

B^Yl ^'""^^Bj , Bj G S(m) . (2.11) 

j=0 

In all of the above symbol classes the composition of the corresponding Weyl operators is 
well defined in the following sense (see, e.g., [pS99|| ): 



Lemma 2.2. Let mi,m2 be order functions. For Bj G S(mj) the product of the associated 
Weyl operators reads in terms of their symbols 

B,B2 = B^2, (2.12) 

where {Bi,B2) ^ Bi^B2 is a bilinear continuous map from S(mi) x S(m2) to S(mim2). 
It is explicitly given by 

{B,i^B2){p,x) = et-(^-^^^^-^^)5i(p,a;)i?2(g,t/) , (2.13) 

y=x 

where cr{vp, v^', Wp, w^) '■= ■ Wp — Vp ■ Wx denotes the symplectic two-form on x M.'^. 
Furthermore, the asymptotic expansion of ( ^.iSj ) in S(mim2) reads 



^i#52)(p, ^) ~ XI fc! ( ^^(^p' ^f) ) ' ^) ^2(q, y 

k=0 ' ^ ^ 



(2.14) 



q=p 
y=x 



The Pauh Hamiltonians ( |2.2| ) that we are going to study below can be viewed as Weyl 
quantisations of hermitian symbols H G S°i(m), where m is an order function with m > 1. 
The symbol H then has an asymptotic expansion of the type ( 2.11| ), 

oo 

H WHj , Hj G S(m) , (2.15) 

and the principal symbol Hq is supposed to be a scalar multiple of !„, i.e. Hq = ifo.sln. 
Since in quantum mechanics Hamilton operators must be self-adjoint and bounded from 
below, we require the real valued function ifo,s to fulfill the following properties: 

1. Hq s is bounded from below. 

2. There exists an energy value E eM. and some e > such that HqI{[E — e, E + e]) C 

X is compact. 

3. Hq s has no critical value in [E — e, E + e]. 

4. {Hq s + i) is elliptic in the sense that 

\Ho^s{p,x) + i\>Cm{p,x), (2.16) 

with some C > 0. 

If h is small enough, the above properties of iJo,s ensure that the spectrum of H is purely 
discrete in any interval which is properly contained in [E—e, E+e], and that H is essentially 
self-adjoint on C^(M'^) ® C". For simplicity, we denote the self-adjoint extensions, which 
are the relevant quantum mechanical Hamiltonians, also by H. Furthermore, if / G C^(M), 
the operator f{H) defined by the functional calculus given by the spectral theorem is a 
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Weyl operator with symbol G S°(m ^) for every r G M. Its asymptotic expansion is given 
by 



oo 



~Z^^V. , (2.17) 

j=0 



with 0o(p, a?) = f{ Ho sip,x))ln a nd 0i(p,a;) = Hi{p,x)f'{Ho^s{p,x)). The above prop- 
erties are proven in |[Rob87| , PS99|| for the case of scalar valued symbols. These proofs can 



immediately be carried over to the present situation in which the principal symbol Hq of 
the Hamiltonian is a scalar multiple of 1„. 

We remark that the requirement Hq G S(m) together with (|2.16| ) implies that 



C'm\p,x) < 1+Hl^{p,x) < D'm\p,x) + 1 , (2.18) 

so that (l+i^Q g)^/^ is an order function. Moreover, the condition Hj G S(m) is equivalent to 
Hj G S((l + iJQ j,)^/^), which is the requirement on the symbol imposed in [[Rob87| , [HMR87|| . 



Below we will study the semiclassical behaviour of quantum mechanical observables 
and of their expectation values in eigenstates of the Hamiltonian. For these purposes it 
will be advantageous to restrict attention to bounded observables. In order to characterise 
a sufficiently large class of such operators we employ 

Proposition 2.3 (Calderon-Vaillancourt). Let B G S°(l), then the Weyl quantised 
operator B is hounded on L'^iW^) ® C". Moreover, for h G (0,1] there exists an upper 
bound for the operator norms of B. 

The original version of this Proposition goes back to Calderon and Vaillancourt [ pV71|| . 



In the form presented here, it can be found (for scalar valued symbols) in |PS99|| ; that proof 
can be directly carried over to the present situation. Furthermore, if we require the symbols 
G S°(l) to be hermitian matrices, their Weyl quantisations B are symmetric operators 
that can be extended to self-adjoint operators on all of L'^(W^) ® C"; the latter we also 
denote by B. The above properties can immediately be generalised to the case B G S'^(l). 
Then, however, one has to extract a factor of from the norm of B in order to obtain 
the bound on the norms when h varies in (0, 1]. 

3 Semiclassical time evolution and Szego limit for- 
mula 



The quantum mechanical time evolution as, e.g., governed by the Pauli equation ( |2.1| ) 
requires to investigate the Cauchy problem for the operator 

ifi|-// (3,1) 
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on L'^{M.'^) ® C". This problem can be solved by introducing the one-parameter group of 
unitary operators 

U{t) := e"^^* , (3.2) 

which are well defined for alH G M since H is (essentially) self-adjoint. The time evolution 
of a quantum mechanical observable B, which arises as the Weyl quantisation of a symbol 
B G S^i(l), is then given by 

Bit) -.= 11^1)3 U it) . (3.3) 
The time-evolved observable ( p.3|) therefore satisfies the Heisenberg equation of motion 

|5(t) = ^[i^,5(t)] , m = B. (3.4) 

For Weyl operators H and B with scalar symbols it is well known that for finite times t the 
propagators U (t) are semiclassical Fourier integral operators which arise as quantisations of 
the Hamiltonian fiow generated by the principal symbol Hq s of the quantum Hamiltonian 
H . In addition, the time evolution ( |3.3| ) respects operator classes in the sense that if 13 is 
a quantisation of a symbol in some suitable class such as S^i(l), B{t) is again an operator 
with symbol in the same class (see, e.g., |[Rob87| |). 



In the case of matrix valued symbols the situation is different; in physical terms this has 
to do with the need to propagate the internal (i.e. spin) degrees of freedom in addition to 
the translational ones. To begin with, let us hence introduce the Hamiltonian vector field 
Xho^s •= {—dxHo,s, dpHo^s) associated with the scalar factor ifo,s of the principal symbol 
Hq, and denote by $* the fiow generated by ^Hq.s- That is, {p{t),x{t)) = $*(p, cc) is a 
solution of Hamilton's equations of motion 

Pit) = -d,Ho^s{p{t), xit)) , xit) = dpHo^sipit), xit)) , (3.5) 

with initial condition (p(0),a;(0)) = {p,x). This describes the classical dynamics of the 
translational degrees of freedom, whereas the propagation of the spin degrees of freedom 
is only contained in the dynamics of the (matrix valued) symbol B{t). For the following 
we suppose that i? is a Weyl operator with symbol B G S^i(l). Then ( |3.3|) yields the 
observable B [t) , which is a Weyl quantisation of some symbol B it) that in general will not 
be in the class S^i(l), although, however. Bit) clearly remains a bounded operator. We 
hence base the following construction on the formal asymptotic expansion 

oo 

i?(t)~5^r^+'i?Kt) , (3.6) 

1=0 

whose coefficients can be determined from the Heisenberg equations of motion ( p.4| ), once 
these have been transfered to the level of symbols with the help of the product formula 



8 



( p.l4| ); the latter also applies to operators with symbols that have a formal asymptotic 
expansion of the type ( p.6|) . We thus obtain the recursive Cauchy problem, 



d_ 
di 



Bi{t)-{Ho,Bi{t)}-i[H,,Bi{t)] 



0<k<l~l 
j+\a\ + \l3\=l-k 



a 



with -B;(0) = Bi, compare |[lvr98| , ch.2.3]. Here 

{A B} := dpA ■ d^B - d^A ■ dpB 



(3.7) 



(3J 



denotes the Poisson bracket of A,B e C°°{R'^ x M'^) (g)C"^"; notice that in general {A, B} ^ 
— In leading semiclassical order (p.7|) now yields the following equation for the 
principal symbol Boit), 



^Eo(t) - {Ho, Boit)} -i[H,, Boit)] = , 
which for Hamiltonians with a scalar sub-principal symbol Hi reduces to 

|i?o(t) - {Ho, Bo{t)] = A (Boit) o = . 



In order to solve the general case, (|3.9|) is rewritten as 
d 



dt 



[d-\p, X, -t) Bo{t){^-\p, x)) d{p, X, -t)] = 



where d has to fulfill 

d{p,x,t) + iHi{^\p,x))d{p,x,t) = 



d{p, X, 0) = 1„ 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



Here the time derivative has to be understood along the trajectory $*(p, x). The quantity 
d already appeared in | [BK99a| , [BK99b| j, where it was introduced in order to describe the 
(semiclassical) propagation of the spin degrees of freedom along the trajectories of the 
Hamiltonian flow $*. 



We are now in a position to calculate the principal symbol of B(t) from ( 3.11 ), 

Bo{t){p,x) = d{<!>\p,x),-t)Bo{<!>'{p,x))d-\<!>\p,x),-t) . (3.13) 
Equivalently, employing the property 

d{^\p,x),-t) = d-\p,x,t) (3.14) 
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that can be deduced from ( p.l2| ) (see, e.g., ||BN99|| ) one obtains 

Bo{t){p,x) = d-\p,x,t) Bo{<^\p,x))d{p,x,t) . 



(3.15) 



Notice that in ||BM 99|| the quantity r{t,p,x) = d^^{p, x,t) is used instead of d{p,x,t). 
If we take into account that d{p, x, t) is unitary, which follows from ( p.l2| ) since Hi is 
hermitian (see, e.g., PK994 PN99|| ), then (|3.15|) also reads 



Bo{t){p, x) = d\p, X, t) Bo{^\p, x)) d{p, X, t) 



(3.16) 



which is the form that we will use below. 

In principle one could determine all coefficients Bi{t) in this fashion, but one must be 
aware of the fact that the resulting symbol B{t) will in general not be in S^i(l), unless 
one restricts the growth of the symbol H in such a way that one would exclude the Pauli 
Hamiltonians ( |2.2| ) that we are interested in. If, however, one restricts attention to ob- 
servables whose symbols have compact support, B G C^(M°' x R"') C"'^"^, this problem 



is avoided and one obtains a Weyl symbol Bsum{t) G C, 



oo / 
I 



® C"^" by, say, Borel 



summation of the asymptotic expansion ( p.6|) , with the solutions Bi{t) of the recursive 



Cauchy problem (|3.7|) entering. On the level of the corresponding operators one therefore 
finds 



< Ch 



N 



for all AT G N 



(3.17) 



where || ■ ||£ denotes the norm on C{L'^{W^) ® C", L'^{W^) ® C"). 

The same result, which is a variant of the Egorov Theorem [ [Ego69|| , can be achieved 
for general B G S^j(l) under suitable requirements on the symbol of the Hamiltonian H 
|v?9|, ch.2.3]): 



see 



Proposition 3.1. Let H he in S(m) and B G S^i(l) and suppose that 



\d"dl^Hj{p, x)\\<C for all (p, x) G M"' 



X 



and |a| + +2j > 2 . (3.18) 



Then the estimate 



Bsumit) - e^""' B e 



<Dh 



N 



(3.19) 



holds for arbitrary N > and t G [0, T] . 



Let us remark that under the assumption (|3.18|) the Hamiltonian vector field X}j^ ^ 
grows at most linearly at infinity. Therefore, a trajectory $*(p, x) cannot blow up at finite 
times so that the fiow exists globally on x (see, e.g., ||Rob87|| ). We have not made any 
attempt at improving the bounds on the time T, in terms of up to which Proposition 
3.1 holds, since for our further purposes we are mainly interested in the relation ( 3.16| ). 
For the scalar case such improvements have, however, been established in | |BGP99| , [BK99 |. 

As a second essential input for the proof of quantum ergodicity, in addition to the 



Egorov property ( |3.16D , we require a Szego limit formula (see, e.g., ||Gui79| |) that connects 
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averaged expectation values of an observable semiclassically with a classical average. On 
the quantum mechanical side we consider an interval I{E,h) := [E — huj,E + hw], with 
some a; > 0, such that I{E, K) C [E — e,E + e] ii h is sufficiently small. According to 
the assumptions 1.-4. on the Hamiltonian H in section ^ the spectrum of H in I{E,h) 
is therefore discrete. We then denote the number of eigenvalues contained in I{E, K) by 
Ni, and {ipk} shall be the (orthonormal) eigenvectors of H associated with the eigenvalues 
Ek G I{E, h). On the classical side, let VLe '■= HqI{E) be the hypersurface of energy E in 
phase space and denote the normalised Liouville measure on Qe by d/x^, i.e. 

dfiE{p,x) = ^—6{Ho,s{p,x) -E) dpdx . (3.20) 
vol \Ie 

In the following we will abbreviate averages of (smooth) matrix valued functions B G 

(:7°°(M^ X R"^) ® C"^" over ^e by 

^e{B):= f B{p,x) dfiE{p,x) . (3.21) 

From now on we also suppose that the Hamiltonian flow $* generated by Hq s is ergodic 
with respect to fiE, i-e. for / G L^{flE,dfiE) and /i£;-almost all {p,x) G Qe 

lim i r f{<i>\p,x)) dt= [ f{q,y) dfiE{q,y) • (3.22) 

This in particular implies that the set of periodic points of $* with periods T > has 
Liouville measure zero. 

The rest of this section will now be devoted to the proof of 

Proposition 3.2 (Szego limit formula). Let H be a quantum Hamiltonian with symbol 
H G S°[(m), where m > 1, that fulfills the properties 1.-4- of section |^. // then B is 
an observable with symbol B G S°[(l) and principal symbol Bq, the following Szego limit 
formula holds, 

lim V {^k,B^k) = - tifiE{Bo) . (3.23) 

Proof. In principle we adopt the method presented in | pS99| | , and modify it appropriately 
where necessary. Let us ffist recall that if g E C^(M) is suitably chosen, with g{X) = A on a 
neighbourhood of the interval I{E, h), the operator g{H) has the same spectrum in I{E, h) 
as H itself. Furthermore, the symbol Hg G S° (1) of g{H) has an asymptotic expansion 
that coincides on HqI(I(E, h)) with that of H. Since below we are localising in energy to 

the interval I{E, H) so that we can consider g{H) instead of H, from now on we simply 
suppose that H G S° (1). 
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Let now x ^ C'^(]R) be given with x = 1 on I{E,h) and such that the spectrum of 
H in supp X is discrete. This might require to choose h small enough. Then consider the 
(energy localised) quantum mechanical time evolution operator 



U^{t) :=e~T.^'x{H) . (3.24) 

Up to an error of order 0{h°°) in trace norm, this operator can be approximated by a 
semiclassical Fourier integral operator with kernel 

K^{x, t) = -i— / an{x, y, t, ^) et(5(-'«'*)-«-^) d^ , (3.25) 

if \t\ is small enough (see, e.g., |PS99|| ). We remark that here the amplitude an takes values 
in C"^", thus it also represents the internal (i.e. spin) degrees of freedom. As explained 
BK99a|| , the phase S in (p.25|) is then given as the solution of the Hamilton- Jacobi 



equation 

Ho,s{d..S{x,tt),x) +dtS{x,tt) = , S{x,tO)=x-^. (3.26) 

In leading semiclassical order the transport equation for aj^ ~ + hai + . . . can be 
solved by a separation of the internal degrees of freedom from the translational ones. 
The latter lead to the expression known from the scalar case (see, e.g., |PS99|| ), whereas 
the modifications required by the matrix character of Oq are provided by the solution 
d{p,x,t) of (|3.12| ), see ||BK99a|| . For the present purpose, however, one only needs the 
initial condition afi\t=o = x(-f^o,s)ln + 0{h). 

In a next step we consider p G C°°(M) with Fourier transform p G C^(M) such that 

Tr ^ ^ m e^^* B U^{t) = Y. ^(^'^) ^^'^) P {^^Y^) ' ^^.27) 

where Tr(-) denotes the operator trace on L'^{M.'^) C". Approximating U^{t) with the 
help of ( |3.25|) , in leading semiclassical order one then has to calculate 

— i— / / / p(t)tr(i?o(9.^,a;)ao(a.,x,t,0)e^(^(^-«'*)-«-^+^*)deda;dt (3.28) 

with the method of stationary phase. The stationary points of the phase S{x, ^, t)—$,-x-\-Et 
are given by (^^j, Xgt, tgt) G M'^ x M'^ x R such that (^g^, Xgt) G Qe is a periodic point of the 
Hamiltonian flow $* with period tgt- Since we have assumed that E is not a critical value 
for Hq s, the periods of on Qe do not accumulate at zero, see |[Rob871| . Hence, if the 
support of p is chosen small enough, the manifold of critical points contributing to ( ^.28 ) 
is given by fi^; x {0}. In analogy to ||BK99a|| we thus obtain 

J2xiE,){^,,m)pi^^^^=x{E)~^^^{^^ . (3.29) 
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Moreover, since we require ifo,s to be such that no E' G [E — e,E + e\ is a critical value and 
all VIe' are compact, ( p.29|) holds true with E replaced by E' uniformly in [E—e, E+e\. In a 



last step we now apply the Tauberian Lemma of ||BPU95|] , which takes non- negative weights 
into account. To this end we have to ensure that {ipk, Bipk) is non-negative. However, 
since B is bounded, this can always be achieved by adding a suitable constant. Moreover, 
according to our above choice of x we find that x{B) = 1 and x{Bk) = 1 for all E^ G I{E, h). 
Therefore 

^ (^,,5^,) = ^^^tr/xM5o) + o(ni-'^). (3.30) 

Repeating the above reasoning with the identity instead of 5, one can express the number 
Nj of eigenvalues in I{E, K) semiclassically, 

Ni = - + o(/^-^) . (3.31) 



Thus, ( |3.3CI| ) and ( |3.31| ) together finally yield the Szego limit formula (|3.23| ). □ 



4 Quantum ergodicity 

In the following we will restrict attention to the case of spin s = 1/2, which is both 
the simplest and physically most important situation. This means that below n = 2s + 
1 = 2 will be chosen, so that all symbols of observables and Hamiltonians take values in 
the hermitian 2x2 matrices. We also assume that the quantum Hamiltonian describes 
the coupling of translational and spin degrees of freedom as in ( |2.2| ). This restricts the 
subprincipal symbol Hi of to be a traceless hermitian matrix. As a consequence, the 
spin transport equation ( ^.12] ) is solved by a spin transport matrix d{p,x,t) G SU(2). 

Our strategy of approaching quantum ergodicity is inspired by the method introduced 
in [ Zel96 , ZZ96|] , which does not require to rely on a positive quantisation, such as anti-Wick 



or Friedrichs quantisation. It is rather based on an analysis of the expression 

S2{E,h):=^ J2 \{4^k,B^k)-ltTfiE{Bo)f , (4.1) 
^ Ei,ei{E,h) 

which is the variance of the expectation values of the quantum observable B about the 
classical mean value of its principal symbol Bq. We are in particular interested in the 
behaviour of ([4.1|) in the limit h 0. For this purpose we introduce the bounded and 
self-adjoint auxiliary operator 



BT:=^j^ U\t)BUit) dt-^tifXE{Bo 



(4.2) 
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Its expectation values in eigenstates of the Hamiltonian are 

{ipk,BTtpk) = {ipk,BtlJk) - ^tTfiE{Bo) , (4.3) 

such that 

S2{E,h) = ^ J2 \{i^k,BTi^k)f . (4.4) 

Using the Egorov property (|3.16|) , the principal symbol of the auxiliary operator reads 

1 1 
BT,oip,x) = -J d\p, x,t) Bo{^\p, x)) d{p, X, t) dt- -trnEiBo)l2 . (4.5) 

Comparing this expression with the analogous one obtained in the scalar case, one observes 
that the principal symbol Bq is not only transported by the flow $*, but also conjugated 
with the spin transport matrix d{p, x, t). Thus both the translational and the spin dynam- 
ics are involved. This observation suggests that the ergodicity of the flow $* will no longer 
suffice to yield quantum ergodicity. One therefore has to combine the classical dynamics 
of the translational degrees of freedom and the spin dynamics in a suitable way, and then 
one demands ergodic properties of the combined dynamics. In order to achieve this we 
employ the same construction as in | [BK99b| | and hence introduce the product phase space 

A< := fiij X SU(2) . (4.6) 

The combined flow : ^ is then defined as an SU(2)-extension of the Hamiltonian 
fiow $* on VIe, i.e. for (p, x) G VIe and g E SU(2) we set 

y*((p, x),g) := ($*(p, x), d{p, X, t)g) . (4.7) 

The initial condition = id is obviously fulfilled, and = o follows from the 
composition law 

d{p, x,t + s) = rf($*(p, x), s) d{p, X, t) (4.8) 

that derives from the spin transport equation ( p.l2| ). In ergodic theory such a combined 
dynamics is also known as a skew product (see, e.g., ||CFS82|| ). On M. the product measure 
H := fiE >^ f^H consisting of the Liouville measure fiE on Qe and the normalised Haar 
measure fin on SU(2) is introduced. This measure is normalised and invariant under the 
fiow y*, since /^s is normalised and invariant under and hh is both left and right 
invariant. Ergodicity of Y^ on M. with respect to /x then means that for F e L^(A^,d/i) 

hm i r F{Y\{p,x),g)) dt = [ F{{q,y),h) dfi{{q,y),h) (4.9) 

holds for yU-almost all {{p,x),g) G Ai. In particular, if one chooses F to be independent 



of g, then (|4.9|) reduces to the condition (|3.22|) , so that the ergodicity of the extended fiow 
y* on Ai implies the ergodicity of the fiow $* on the base manifold Qe- 

Our main result now states the effect of the ergodicity of Y^ on the semiclassical asymp- 
totics of eigenfunctions of H. 
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Theorem 4.1 (Quantum ergodicity). Let H be a quantum Hamiltonian with symbol 
H G S°i(m), where m>l, and principal symbol Hq = i/o,sl2; which fulfills the conditions 
1.-4- of section^ and, furthermore, satisfies 

\\d'^d^Hj{p, < C for all {p, x) eM.'^ x M.'^ and \a\ + \l3\ + 2j > 2. ( glSD 

Then, under the condition that is ergodic on M. with respect to fi, in every sequence 
{ipk I Ek G I{E,h)} of orthonormal eigenf unctions of H there exists a subsequence {ipkj \ 
Ek^ G I{E,h)} of density one, i.e. 

,. *U \ Et, e I{E,h)] 

te5#{*|B,e/(B,s)}=l- (^-l"' 

such that for every quantum observable B with hermitian symbol B G S[!j(l) and principal 
symbol Bq 

lim {i^k.^BiJk,) = ItifiEiBo) . (4.11) 

Moreover, the subsequence {^kj] can be chosen independent of the observable B. 

Proof. An application of the Cauchy-Schwartz inequality on the right-hand side of (|4.4| ) 
yields an upper bound for the quantity ( ^?T| ) that reads 

S2{E,h)<^ J2 i^kABTf^k). (4.12) 

^ Ek£l(E,h) 

In order to determine the limit as h ^ oi ( |4.12D one can now apply Proposition p.2| , 
which gives 

hmS2{E,h) < ltTfiE{iBT,o?) . (4.13) 

Quantum ergodicity then follows, if the bound on the right-hand side can be shown to 
vanish. This will indeed be possible in the limit T — oo. In order to achieve this we first 
employ the ergodicity of in that we choose F{{q, y),h) = Bq^q, y) h in the relation 
( ri) , and thus find 

lim / g^dAp,x,t)Boi<^\p,x))d{p,x,t)g dt= [ fiE{Bo) h dfinih) (4.14) 

T^oo I Jq JSU(2) 

for /i£;-almost all {p,x) E Qe and /ij^-almost all g G SU(2). In terms of the principal 
symbol ( [4.5|) of the auxiliary operator Bt this means 

lim g^ BT,Qip,x)g = I fiE{Bo) h d^inih) - ^tr fiE{Bo) I2 ■ (4.15) 
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Our next goal is to calculate the right-hand side of ( [4.151 ). To this end we represent the 
hermitian 2x2 matrix ^e{Bq) as a linear combination of I2 and the Pauli matrices 0"^, i.e. 
fiE{Bo) = ^tT fiE{Bo) + b ■ (T, where b e M^. We then recall that for every g e SU(2) the 
adjoint map Adg : su(2) su(2) is defined as Adg{X) = X g. Thus, Ad/, (6 ■ a) e su(2) 
can be expanded in terms of the Pauli matrices ak such that b ■ a h = [ip{h)b) ■ cr. The 
map : SU(2) S0(3) that results in this way can be identified as the universal (twofold) 
covering of SO (3) by SU(2). Therefore 



/ flE{Bo)hdflH{h) = ltTfiEiBo)l2+ ( [ ^{h)dflHih)b 

JSV{2) ^ \JSV{2) 



(4.16) 



We now show that the second term on the right-hand side of ( |4.16|) vanishes. To this 
end we multiply the integral over SU(2), which yields a 3 x 3 matrix, with an arbitrary 
orthogonal matrix R G SO (3) from the left. We then exploit the fact that there exists 
g G SU(2) such that R = ip{g), together with the left invariance of the Haar measure, in 
order to conclude that 



R I ip{h) dfinih) = / ^pigh) duuih) = / (p{h) d^nih) . 

'SU(2) JSU(2) JSU(2) 



(4.17) 



Since hence the 3x3 matrix represented by the integral over SU(2) is invariant under left 
multiplication by an arbitrary element of SO (3), it must be the zero matrix. Therefore, 
the right-hand side of ([4.15|) vanishes. 

We now choose some g G SU(2) such that ( [4. 141 ) holds and therefore obtain that 



lim tr(5r,o(p, x)f = lim ii[g'^ Br^iP, x) g^ = 



(4.18) 



holds for /i£;-almost all {p,x) G Qe- Thus, after an integration over Qe, this together with 
( [4.131 ) implies that 



\imS2iE, h) = . 

ft— >o 



(4.19) 



According to a standard argument in the proof of quantum ergodicity for scalar Hamil- 
tonians, see ||Zel87[ , |CdV85[] , the vanishing of S2{E,h) in the semiclassical limit implies 
the existence of a density-one subsequence {i'kj}jm C {i'k}k<=N such that ( [4.11[ ) holds. 
Another standard, diagonal construction then ensures that a subsequence can be chosen 
that is independent of the observable, see ||Zel87[, [CdV85[|. □ 



5 Discussion 

In the case of scalar Hamiltonians quantum ergodicity is often interpreted in terms of 
Wigner- or Husimitransforms of eigenfunctions. In this context one concludes that along a 
subsequence of density one the Wigner- or Husimitransforms of eigenfunctions of a quan- 
tum ergodic Hamiltonian weakly converge, as distributions or measures, respectively, to 
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Liouville measure. Thus the hfts of eigenfunctions to phase space become equidistributed 
on the hypersurface Qe of energy E. In the present situation an analogous interpretation 
first requires to introduce matrix valued Wigner- and Husimitransforms. In general, the 
Wignertransform of ip E S'iM.'^) ® C" is given by 

f i — 1 1 
W[^]ip,x)= e-T^P-y^[x--y)(^i){x + -y) dy . (5.1) 

Then expectation values of observables B with symbols B G S'^(l) in states described by 
i\) e L2(R'^) ® C" read 

{^,Bij) = -1- / / tT{wmp,x)B{p,x)) dpdx . (5.2) 



In order to convert the statement of Theorem 4.1 into one about the matrix components of 



Wignertransforms we introduce the special observables with symbols = bErs G 
S°(l), where 6 is a real valued function on phase space that is independent of h, and the 
four constant matrices Ers are defined by 

^"=(o o) ' ^''=[0 o) ' ^'' = [1 o) ' ^''={o i) • ^^-2) 



Although the off-diagonal symbols b E12 and b E21 are non-hermit ian, Theorem ^]T| can be 
applied to all of the observables B^^^^ in an obvious manner. Together with the relation 
( ^.2|) this then reveals that along the subsequence {i'kj} of density one specified in the 
Theorem the matrix components of the Wignertransforms of the eigenfunctions weakly 
converge, as distributions on C^(R'^ x R"^), to either Liouville measure or to zero. More 
specifically. 

Urn —^W[iJk] = -^7^S(Ho,- E)l2 (5.4) 

component- wise in ViW^ x M*^). In particular, the semiclassical limit of the Wignerdistri- 
butions is a scalar multiple of the identity matrix. This means that the 'spin up' and 'spin 
down' components become identical and equidistributed over the hypersurface of energy 
E in phase space. Moreover, there occurs no mixture between 'spin up' and 'spin down' 
components, which can also be seen on the right-hand side of (|4.11| ) since there only the 
diagonal elements of the principal symbol Bq contribute. 

Upon introducing matrix valued Husimitransforms through 

/ / Vr[V^](q,y)e-^«P~'')'+("-y)') dgdy , (5.5) 



Hmp,x) 



one can also consider anti-Wick quantisations of symbols B G S°(l). Expectation values 
of the corresponding anti-Wick operators Baw then read in analogy to ( p.2| ) 



Oj^Baw^) = j7^, [ [ tT{H[,p]{p,x)Bip,x)) dpdx . (5.6) 
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Since, moreover, for B G S''(l) 



B -B 



AW 



0{h) , (5.7) 



one can replace B by B^w in ( [1.11| ) so that Theorem |4.1j imphes an analogue of (|5.4| ), i.e. 



li"^ TT^^t^'^J dpdx = J I2 d/XE , (5.8) 

where the convergence has to be understood component-wise as a weak convergence of 
probability measures on the phase space M.'^ x M.'^. 

As a further point we now want to discuss whether it is actually necessary to introduce 
the SU(2)-extension of the Hamiltonian flow and to require ergodicity of in order 
to obtain quantum ergodicity. To this end we provide an example showing that in general 
it would not suffice to demand only ergodicity of Let us therefore introduce a quantum 
Hamiltonian H of the type introduced in sections 0. - ^. with symbolQ 

H{p, x) = Ho,sip, x) I2 + hC{p, x) aj , (5.9) 

where aj is one of the Pauli matrices and Hq s shall be chosen such that is ergodic ohQe- 
Now cTj can obviously also be considered as a bounded self-adjoint operator on L^(]R'^) (8)C^ 
which commutes with H; in fact, ^aj is the j-th component of the spin observable. Hence, 
one can introduce joint eigenvectors ipk G L'^(M.'^) (g) of and aj. Since = I2, the 
eigenvalues of aj are Afc = ±1. Introducing Uj G U(2) such that UjCTjUj is diagonal, we 
can switch to eigenvectors ipk := Ujipk such that 

= (^^''^^^ if Xk = +1 and = (^J-)^ if = -1 . (5.10) 

The expectation values of an observable B in these eigenvectors (pk therefore read 



(y.l+\5n^l+)) ifAfc = +l, 
B,,^[-^) if A, 



{^,,B^,) = {)^t_;-''^f_,i (5.11) 



where (■, •) denotes the scalar product in L'^(W^). Since the subsequences of the ipkS with 
Xk = +1 and Afc = —1, respectively, are each of density one half, quantum ergodicity 
cannot hold for general observables B. 

On the other hand, the above example does not fulfill the requirements of Theorem ^?l| 
so that no contradiction occurs. To see this consider the spin transport equation ( |3.12| ), 
which in the present case reads 

d{p,x,t) + iC{<^>\p,x))ajd{p,x,t) = , d{p,x,0) = l2, (5.12) 



-'^We owe this example to Stefan Keppeler. 
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and is solved by the expression 



d{p, X, t) = cos{a{p, X, t)) I2 — i sin(Q;(p, x, t)) crj , (5.13) 

where 

a{p,x,t)= [ C{<^'{p,x)) ds . (5.14) 
Jo 



The exphcit solution ( |5.13D of the spin transport equation demonstrates that, even if the 



flow $* on the base manifold fi^; is ergodic, its SU(2)-extension cannot be ergodic on the 
product phase space A^. This is due to the fact that with ( ^.13|) one only explores a one 



dimensional submanifold of the three dimensional group manifold of SU(2). We therefore 
conclude that in the case of Pauli Hamiltonians ergodicity of alone is not a sufficient 
criterion for quantum ergodicity. 
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